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Abstract 

We consider leptogenesis in scenarios with many neutrino singlets. We find that the 
lower bound for the reheating temperature can be significantly relaxed with respect 
to the hierarchical three neutrino case. We further argue that the upper bound for 
the neutrino mass scale from leptogenesis gets significantly lifted in these scenarios. 
As a specific realization, we then discuss an extra-dimensional model, where the large 
number of neutrinos is provided by Kaluza-Klein excitations. 
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1 Introduction 



Among the many open issues within astro-particle physics and cosmology, the origin of the 
observed baryon asymmetry of the universe (BAU) [1] 

77^^^= 6.ll°i X 10-^° (1) 

Th-y 

is one of the most frequently discussed, as it cannot be explained within the standard 
model of particle physics (SM). Here, = nf, — is the difference of today's baryon and 
anti-baryon densities and is the corresponding value for the photons. 

Especially, in the context of inflation, which can solve many cosmological problems, such 
as the flatness and the horizon problem, one typically starts with an empty universe at the 
end of this epoch, which implies that somehow a baryon- asymmetry needs to be produced 
before big bang nucleosynthesis (BBN) sets in. 

One of the most popular models that can produce this asymmetry is leptogenesis [2-9] . The 
popularity of this model can be ascribed to the fact that it does not need any additions to 
the SM other than right-handed neutrinos, which are also needed to make the SM neutrinos 
massive, as indicated by their observed oscillations. If these right-handed neutrino states 
are very heavy, they can nicely explain the smallness of the observed neutrino masses via 
the sec-saw mechanism. In leptogenesis, these state are indeed assumed to be heavy. After 
their excitations in the early universe they decay in a CP violating manner, leaving a 
lepton asymmetry, which is then partially transformed to a baryon asymmetry by so called 
sphaleron processes [10, 11]. Lately, much attention has also been given to possible flavor 
effects in this scenario [6, 12-15]. 

However, the typical lower bound for the reheating temperature > lO^GeV [8,9,16], which 
is needed to produce a large enough baryon asymmetry via leptogenesis is only marginally 
compatible with the upper bound on the same value in supergravity theories Tr < 10^ — 
lO^^GeV [17, 18], where an exceeding of this temperature would alter the predictions of BBN 
due to an overproduction of gravitinos. Bounds for more constrained models can be found 
in Refs. [19-25]. 

Several ways to overcome this discrepancy have been considered, such as resonant leptoge- 
nesis [12] and soft leptogenesis [26,27]. In the present approach, we consider the effect of 
many heavy right-handed neutrinos (i.e. SM gauge singlets) on standard thermal leptogen- 
esis in the one-flavor approximation. We flnd that a large number of heavy singlet states in 
the weak wash-out regime can signiflcantly lower the required reheating temperature with- 
out any resonant enhancement of the CP asymmetries. For the strong wash-out regime we 
find a slight decrease of the final baryon asymmetry. While the strong wash-out scenario is 
typically considered as more natural in the standard case, we argue that the weak wash-out 
regime can be more natural in models with many singlets. We also consider the limits of this 
scenario, which are reached once the heavy singlets come to dominate the energy density of 
the universe. 

Extensions of the standard model with many singlets can be motivated in several ways. 
Here, we focus on an extra-dimensional model that includes Kaluza-Klein excitations of 



1 



additional singlet states. These excitations can be considered as right-handed neutrinos. 
The model is found to naturally bring all the ingredients that are needed to lower the 
required reheating temperature within the present approach. Other motivations for setups 
with many singlet neutrinos come from string models [28,29]. 

The structure of the paper is as follows. After a review of standard thermal leptogenesis in 
section 2, where we also derive some rule-of-thumb formulae for the final asymmetries and 
present parameter bounds for the scenario, we consider possible effects of many heavy singlet 
states on leptogenesis in section 3. Here, the simplified considerations from the previous 
section will help us to understand the effects of the new particles, while we also present 
numerical results that support our considerations. In section 4, we present the already 
mentioned extra-dimensional scenario that contains all the required ingredients to yield the 
observed BAU at lower reheating temperatures. We also consider bounds for this scenario 
from production and decay of extra-dimensional graviton modes. In section 5, we then 
conclude with a brief summary of our results and a short outlook on possible continuations 
of the present work. In Appendix A, we present the reaction rates used for the numerical 
calculations presented throughout the paper, whereas we provide the reader with technical 
background on the Davidson-Ibarra bound in the presence of many neutrinos in Appendix 



2 Standard Thermal Leptogenesis 

In this section we review the standard thermal leptogenesis case [2-9], with three hierarchical 
heavy neutrinos in the one flavor approximation. 

Therefore, we need to add three gauge singlets to the standard model, which are assumed 
to have Majorana mass terms as well as Yukawa interactions with the SM lepton doublets 



Here, we work in a basis in which the Majorana mass matrix is diagonal, while we assume 
Mi ^ Qjki'P) with (0) being the vacuum expectation value (VEV) of the SM Higgs. 

After electro-weak symmetry breaking, the light (SM) neutrinos also gain an effective Ma- 
jorana mass matrix via the well-known see-saw mechanism, with the three eigenvalues mi, 
m2, and ms. To lowest order in (0)/Mj this can be illustrated by the Feynman graph in 



The addition of the singlets also opens new possibilities for CP-violation, e.g. in the decay 
of the heavy right-handed neutrinos, where an interference of tree-level and loop diagrams 
(cf. fig. 2) can lead to different decay rates for the channels N icff and N i(f). Namely, 
in the case of hierarchical heavy neutrinos (Mi -C M2 -C M3) one finds for the decay of the 
hghtest of them [30] 
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Figure 1: This Feynman graph illustrates the well-known see-saw mechanism, where the 
large masses Mj in the propagator of Nj and the VEV of the SM Higgs field lead to an 
effective Majorana mass matrix for the light neutrinos. 




Figure 2: The relevant decay- diagrams for standard thermal leptogenesis, where the in- 
terference of the first diagram with each of the other two can lead to CP violation and a 
net lepton number production per decay. Due to the different projection operators, the wave 
function correction with charge flow in the opposite direction only contributes at higher 
orders ofMi/M^. 



One can show that there exists an upper bound for this asymmetry [16,31], which is called 
the Davidson-Ibarra bound and which is given by 



\e I < (m m ) - ^^^^^'tn. 



that is maximized by mi = 0. In this case we have Aml^^ = m| and 

3 Mima ^n-i ( \ ( -^i ^ 

with almost no freedom left for the size of ma. 

As pointed out in rcf. [15] the possible CP-violation within AL = 1 scatterings can also be 
relevant for leptogenesis. There, it is also shown that the CP-asymmetry produced in these 
scatterings exactly equals the asymmetry of the corresponding decay processes. 

At the beginning of the scenario, we assume to be in a homogenous and isotropic universe 
(possibly after an inflationary period) dominated by the SM particles, which are all in ther- 
mal equilibrium with temperature T^. The initial abundance of the right-handed neutrinos 
is typically assumed to be zero. 

The abundance of right-handed neutrinos at later times is then determined by the differential 
Boltzmann equation [3,32] 
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where z = Mi/T is effectively a reparametrization of time in a background dominated 

universe. Further, A^at^ is the average abundance of the hghtest right-handed neutrino in 
a comoving volume containing one photon at early times, is the corresponding value 
for the equilibrium distribution, Td is their decay rate, and rAL=i the rate for the 2^2 
scattering processes, which violate lepton number by one (e.g. Nii <-> tq). The explicit 
values for these quantities can be found in appendix A. 

The second important equation in this context is the one for the B — L-asymmetry Nb-l 
[3, 15, 32] in the same volume 

= H^D + rAL=i) {Nn, - iV^^J + TwNb-l] , (7) 

where is the reaction rate for wash-out processes consisting of inverse decays and 
AL = 1,2 scattering. The exphcit formulae of the wash-out rates can also be found in 
appendix A. 

A fraction of the final B — L asymmetry will be made up of baryons due to sphaleron 

processes that partially convert a lepton asymmetry into a baryon asymmetry {Nb = 
28/79Nb-l [33,34]). To find the final baryon asymmetry per photon rjB in a late universe, 
one also has to consider the late time production of photons, which leads to a dilution of 
the asymmetry [9]: 

riB ~ 10-^NB-L{t ^ oo) . (8) 

Extensive numerical and analytical analysis of equations (6) and (7) have been performed 
in refs. [2-9, 15]. In this part of the paper, we aim at understanding the behavior of the 
solutions of these equations qualitatively to conceive the effect of many decaying neutrinos 
in the following sections. For reasons of simplicity we will adopt the treatment of ref. [9] 
and restrict our considerations to decays, inverse decays and AL = 1 scatterings involving 
the top quark. We will comment on AL = 2 scatterings and scatterings involving gauge 
bosons later in this section. 

When analyzing the Boltzmann equations of leptogenesis, it is useful to define the parameter 

H{z=l) ' 

which gives a measure for the decoupling of the right-handed neutrinos at the time when 
their equilibrium abundances start being suppressed. We distinguish the weak wash-out 
regime {K ^ 1) and the strong wash-out regime [K ^ 1) and will discuss both of them in 
the following subsections. 

However, before going into detail let us note that the parameter K is often characterized 
by the relation 

K^^, (10) 
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where 

Ml ^ ' 



is the so-called effective neutrino mass and 

Mpi 



10-^eV (12) 



is the called the equilibrium neutrino mass, with the Planck mass Mpi = and G 

being Newton's constant. 

2.1 The weak wash-out scenario 

In the weak wash-out scenario ^ 1) the abundance of right-handed neutrinos is never 
comparable to that of photons and other species due to their weak couplings. For the same 
reason the wash-out processes, which always involve right-handed neutrinos, never come 
into thermal equilibrium (i.e. F^k ^ i/ at all times). 

However, since decays and inverse decays (as well as the corresponding AL = 1 scatterings) 
create B—L asymmetries of opposite sign but equal magnitude, most of the B—L asymmetry 
which gets produced during the decay of the heavy neutrinos only cancels an asymmetry that 
was produced during the population of these states. Therefore, even though the wash-out 
processes are suppressed they are still crucial. 

Let us make some rough estimations for the reaction rates based on dimensional analysis. 
For T > Ml we approximate [32] 

Td ~ ^iD^^{gg%Mi^ (13) 



^-^{99%T%'lTn (14) 



AL=1 



where gtop is the Yukawa coupling constant of the top quark. 

Since the right-handed neutrinos will never reach equilibrium abundance until this value 
becomes strongly suppressed we can put — N^^ for 2; < 1. Using eq.(6) we can 

therefore approximate 

, (rp + Fal=i)a^^, ^ 

NnAz^^)^ fw -^K (15) 

H{z = 1) 

with A?"^ = 3/4 (at early times). 

Without wash-out effects, the B — L asymmetry around this time would be Nb-l{z = 1) pa 
—EiNNiiz ^ 1). However a small fraction of this has been washed-out. To find the size of 
this fraction, we further approximate 

~ l^ID + l^AL=l ~ , (16) 
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Figure 3: A typical example for a weak wash-out scenario. The green/dashed line shows 
the equilibrium density of the right-handed neutrinos Ni per relativistic photon, while the 
red/solid line shows their actual density as calculated from the Boltzmann equations. The 
blue/dotted line shows the absolute amount of the B — L density. The graph was plotted with 
the parameters K = 10~^, Mi = 10^ GeV, Tj. = Mi, ei = 10~^ , and a Higgs boson mass 
of lO^-'^Afi. One can see that the Ni states become mainly populated until z ~ 0.5 and that 
the final B — L asymmetry is much smaller than it is around z ^ 1 due to the discussed 
cancellations. For the final value of \Nb-l\ the numerics give 1.0 ■ lO""*^^, which agrees well 
with our rough estimations in eq.(20). 



where the factor 1/2 is due to fact that half of the inverse decays will finally end up as the 
original states and the factor 2/3 is due to the fact that the s-channel process is suppressed 
more strongly in the weak wash-out before z ^ 1 [3]. 

Using eq.(7) this gives as the very rough estimate (cf. [9]) 

Nb-l{z = 1) ^ eiiV^,(^^l)exp (- ^J"^ (17) 



eiKexp{-K) (18) 
eiK{l-K). (19) 



After the temperature has dropped below Mi, inverse decays and AL = 1 scatterings will be 
suppressed and the particles will decay around H T£,/2, which corresponds to z = K~^^'^. 
The decaying heavy neutrinos will produce an asymmetry, of the same amount as the inverse 
decays and scatterings that populated these states. However, the sign will be the opposite 
and therefore the first term in eq.(19) will be exactly canceled, which shows that we had to 
include the wash-out effects to find a non-zero asymmetry, which then amounts to 

Nb-l{z = oo) ^ . (20) 

A numerically calculated solution can be found in fig. 3, which also illustrates the basic ideas. 
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2.2 The strong wash-out scencirio 



In the strong wash-out scenario {K 1), the wash-out processes come into equihbrium, 
and therefore the B — L asymmetry, which is created during the population of A^i states, is 
washed out before the states decay again. Therefore we do not have a cancellation of the 

two asymmetries as in the weak wash-out scenario. 

On the other hand, the large K also keeps the right-handed neutrino density close to its 
equilibrium value, which means that the majority of the right-handed neutrinos will decay 
while the wash-out processes are still active. This of course reduces the final asymmetry. 

Let use some rough approximations, again. After 2; ~ 1 the wash-out processes will get an 
additional exponential suppression factor, since the right-handed neutrinos now have masses 
larger than the temperature and need to be produced on-shell in all the wash-out processes 
considered in this part of the paper. With regard to eq.(16) we now find (cf. [35]) 

for ^ > 1 which shows that we can neglect the part due to scatterings at late times. 

In a naive approximation the wash-out processes will no longer be in thermal equilibrium 
and therefore not be able to erase a produced asymmetry when ^ 2i7. However, due 
to the rapid suppression of this term, we can make a better approximation. As a result of 
the wash-out processes an asymmetry will behave as 

NB^t) = NB-L(to) exp J' rw(t)d?j = NbM^o) exp ^ ^^dz^ ■ (22) 

As T]^ is very rapidly suppressed, the exponential term can be neglected as soon as the 
integrand is smaller than one. Therefore, a more exact condition for the "time" at which 
the wash-out processes become ineffective is 

^ ^ H (23) 



By this time, the equilibrium number density of right-handed neutrinos has decreased and 
we have 

N{z,) ^ N^'^iz,) ^ ^3/2 exp(-2,)/2 ^ {2z,K)-^ ^ ^ . (25) 

K 

Therefore our final result for the B — L asymmetry is 

iVB_i(z = oo) ^-£1— ^. (26) 

is. 

An example for a numerical solution of the corresponding Boltzmann equations is given in 
fig. 4, which supports our rough estimations. 
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Figure 4: A typical example for a strong wash-out scenario. The green/dashed line shows 
the equilibrium density of the right-handed neutrinos Ni per relativistic photon and almost 
perfectly matches the red/solid line, which illustrates their actual density as calculated from 
the Boltzmann equations. The blue/ dotted line shows the absolute amount of the B — L 
density. This graph was plotted with K = 10^ and all the other parameters as in figure 
3. One can see that the B — L asymmetry from the production of Ni states is almost 
immediately washed-out. For the final value of \Nb-l\ we find 1.5 • lO^""^^, which matches 
our approximative value in eq.(26) roughly to a factor of 3. 



2.3 Upper mass bound for m 



The strong wash-out scenario also places an upper limit on [34,36-38]. Let us illustrate 
the origin of this bound along the discussion in ref. [13], with slightly different numerical 
factors. 

In case of a large m^, all the light neutrinos are degenerate. This also implies [39] 

mi ~ fhi ^ m3 ^ m/ -\/3 and K —- — , (27) 

which tells us that we are in the strong wash-out regime. Thus, the freeze out of AL = 2 
scatterings around z = lis crucial to avoid the wash-out of the generated B — L asymmetry. 
Estimating the rate for these processes to be 



m 



2rp3 



AL=2 



3(0)4 



TV 



{2i 



the condition rAL=2(-2 ~ 1) <^ H{z ~ 1) translates into the upper mass limit 



\ ra 

With eq.(27) we can also transform eq.(4) to 



M < f ^ ) 5 ■ 10^°GeV. (29) 
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Combining this with equations (1), (8), (26), and (27) we find the rough upper mass hmit 
/\/3 < O.leV. In a more precise analysis refs. [9,40] find shghtly higher bounds bounds. 

2.4 Lower mass bound for Mi 

As is well-known, the assumption that the BAU is due to leptogenesis puts a lower mass 
bound on the lightest of the heavy neutrinos (if they are hierarchical). Putting the Davidson- 
Ibarra bound (eq.(5)) together with eq.(8) and our estimates for the B — L asymmetry from 
the Boltzmann equations (eqs.(20) and (26)) we find 

^^A W-'^H^) (loiSv) • wa.h-out; 

1 « ■ "1-"^-' (irSv) (lo^) . ^-"8 --"-out. 

We see that the region around K k, \ will give the largest baryon asymmetry. Reference [7] 
finds 

Ml > 2 • lO^GeV (32) 

as the lowest bound for the lightest of the singlet neutrino masses in case of zero initial 
abundance. 

2.5 Complications 

Several effects have been neglected so far, and we now wish to comment on them: 
AL = 1 scatterings: 

We have neglected AL = 1 scatterings that include gauge bosons. In this paper we adopt 
the treatment from reference [9], assuming that these processes contribute at most in the 
same way as those scatterings that involve the top quark. Being mainly concerned with 
orders of magnitude, this should not change our results. (Also, the neglect of the running 
of the top Yukawa coupling should partially compensate for this.) 

AL = 2 scatterings: 

If these scatterings are in thermal equilibrium at late times, they can wash-out the cre- 
ated asymmetry. Based on dimensional analysis we estimate their rates to be rAL=2 = 
rfP'T^ I (7r^(0)''). The condition that these processes are out of equilibrium around z ~ 1 (cf. 
also [15]) then translates to 

177r3(0)^ 
m^Mpi 

In the weak wash-out scenario these processes can still be in effect around z = \^ but not 
at the time of decay z ~ [2/ KY^"^. The mass bound then changes to 



M < ^ 2 • lO^^GeV. (33) 



Quantum Effects: 

Quantum effects can also play a larger role in leptogenesis scenarios than indicated by our 
treatment. 
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Examples are the running of coupling constants [41] that can affect the final asymmetry 
by an order one factor and thermal quantum effects within the treatment of Boltzmann 
equations [8] that can also have an impact on the final asymmetry. 

Moreover, the semi-classical treatment with Boltzmann equations is of limited validity com- 
pared to the full quantum dynamics described by Kadanoff-Baym equations (see e.g. [42]). 
In spite of recent progress on the subject (e.g. [43,44]) a full quantum mechanical treatment 
of leptogenesis in terms of Kadanoff-Baym equation has not been performed, yet. 

Spectator Processes and Flavor effects: 

The standard treatment for leptogenesis ignores the possible change of chemical potentials 
of particles other than the SM leptons. However, this possible time dependence through 
spectator processes in equilibrium can contribute an order one suppression [45]. 

Lately, much attention has also been given to flavor effects in leptogenesis scenarios, e.g. 
refs. [6,12-15]. The important fact, which had been neglected before, is that at low reheating 
temperatures (T^ < lO^^GeV), the Yukawa couplings for the third (and sometimes second) 
generation leptons are in thermal equilibrium at the time of leptogenesis. This allows us 
to distinguish lepton flavors at that point, and hence one has to use a system of coupled 
Boltzmann equations. Naively, this yields an additional factor of 3 for the flnal asymmetry in 
the weak- and the strong wash-out scenario, while larger discrepancies can also be achieved. 
However, the overall lower bound for Mi cannot be reduced by flavor effects [46]. For 
simplicity, we will only be dealing with the one-flavor approximation in this paper. 

3 Leptogenesis With Many Heavy Neutrinos 

The purpose of this chapter is to examine the effect of additional heavy neutrino states on 
standard thermal leptogenesis especially in the context of the necessary reheating temper- 
ature. The most important result in this context is that this bound can be signiflcantly 
lowered without any resonant enhancement of the CP asymmetries. 

To be more speciflc, we consider a setup with un right-handed neutrinos, of which ries have 
masses that are low enough such that they will reach significant abundances during the 
leptogenesis process (i.e. they have masses not much higher than the reheating temperature). 
In this case, the Lagrangian of this part of the particle spectrum is still given by eq.(2), with 
the difference that the index j is now running from 1 to un, and gf is a tin x 3-dimensional 
matrix. 

The see-saw formula is also not significantly altered, since the diagram in fig.l still gives 
the effective mass of the light neutrinos to leading order (cf. also [47]) 

mu,ik = gJjMr^gjk{(j)f , (35) 
again with the difference that g is now a n^v x 3-dimensional matrix. 

We see that the Yukawa couplings tend to be smaller compared to the standard case with 
similar masses for the singlets, if there are no cancellations between the coupling constants 
or only few dominating couplings. This already implies that the weak wash-out regime can 
be natural in scenarios with many singlets (cf. ref. [29] for a related argument). 
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The Boltzmann equations from the previous section (eqs.(6) and (7)) are easily generahzed 
to 

^ = -;^(rAn + rAL=i,n) (iv^„ - Kl) , (36) 

with 1 < n < rics- One can calculate T^^n and rAL=i,ra as in the standard case, with the 
corresponding interchange of masses and coupling constants (cf. appendix A). 

The second equation now reads 

dNu-L 1 g ^^^^^^^^ ^ ^^^^^^^^ ^^^^ -N^l)+ rw,nNB-L] , (37) 



dz Hz 



n=l 



where can also be calculated as in the previous section, again with the corresponding 
interchange of masses and coupling constants. £„ is discussed in detail in the following 
section and Appendix B. 

3.1 The Davidson-Ibarra bound with many heavy neutrinos 

Let us first examine the CP asymmetry for the decay process of one of the right-handed 
neutrinos. As in the three generation case, we can draw wave-function and vertex corrections 
as in fig. 2. However, once more with the difference that we now have njy of each sort of the 
loop diagrams instead of 3, as we did in the standard case.^ To calculate the CP asymmetry 
we can simply add up the various contributions and find (as in eq.(3)) 

T{N, 4 + $) - T{N, ^ h + $*) 

r{Ni^iL + <^) + T{Ni^iL + <!>*) ^ ' 

_ _ 3 ^ Im[(^^t)2^] M, .^(Mf\ 

3 M 1 . r. + -r. 1 f M, 



Im [{gmW),,] + O j (40) 



167r (0)2 (^^t) . 
with k running from 1 to tiat, now. 

Of course this result is only appropriate if the higher-order corrections in eq.(39) (cf. ref. [48]) 
and resonance effects can be neglected, or if Im [(s'fi'^)^^] vanishes. For simplicity we restrict 
ourselves to this case in the following. 

As we show in Appendix B, eq.(40) leads to the relation^ 

3 Mima .(rn^\(Mi\ 



^In principle, the large number of possible diagrams can invalidate this perturbative approach. However, 
for Yukawa couplings that yield the light neutrino masses without delicate cancellations this is typically not 
the case. 

^The author is very grateful to Yosef Nir for pointing out to him that the Davidson-Ibarra bound does 
lose its validity in the case of many singlets and three light neutrinos, due to the fact that the matrix R 
within the Casas-Ibarra parametrization [49] of the Yukawa couplings is not orthogonal anymore in the 
present case. 
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which matches eq.(5). Yet, in the latter case the CP-asymmetry is maximized for mi = 
and ms is therefore fixed around O.OSeV. In the presence of many singlets is in general a 
free parameter in eq.(41), and therefore the standard Davidson-Ibarra bound is considerably 
weakened. 

However, as we shall see in section 4, the Davidson- Ibarra bound in eq.(5) still applies in 
many cases. This is due to the fact that the number of free parameters for the Yukawa 
couplings can be reduced by symmetries. In particular, the condition of maximal CP asym- 
metries for all decaying particles can motivate such effects. Thus, we will not consider 
possible CP asymmetries beyond the Davidson- Ibarra bound in eq.(5) in the following anal- 
ysis. Nevertheless, it should be kept in mind that there might be further possibilities to 
increase the final baryon asymmetry due to the arguments given in this subsection. 

At this point it is also important to note that the analysis in ref. [15] that determines the 
CP-asymmetry in AL — 1 scatterings is independent of the number of heavy neutrinos. 
Thus, the CP-asymmetries of these processes match the ones of the corresponding decays 
as in the three neutrino case. 

3.2 The strong wash-out scencirio 

Let us first see, how the strong wash-out scenario may be changed by the presence of many 
heavy neutrinos. First, we note that if wc have many neutrinos with ^ 1, it is enough 
to observe the lightest ones, i.e. only those particles that are still abundant in unsuppressed 
numbers, when the wash-out processes due to the hghtest one of the singlets freezes out, 
since all asymmetries produced by heavier particles will have been erased by these processes.^ 
Therefore, we consider the case of ^eff quasi-degenerate neutrinos. 

The addition of extra singlet states has two competing effects on the scenario. The first 
one is the fact that the number of decaying neutrinos is increased, which would increase the 
total amount of the produced B — L asymmetry if it was the only change. However, the 
wash-out rate is also increased, which means the time after which a produced asymmetry 
is not washed-out anymore is later. This effectively reduces the abundance of decaying 
neutrinos. 

Let us make some approximations in the manner of section 2.2. The averaged wash-out rate 
is then (compare eq.(21)) 

- E • l^i99%M, + ^ {gg% t\ .f^e"- , (42) 

with Zi = M/Ti. 

This means the freeze-out "time" of wash-out processes is now defined by the equation 

zl/'e-^-^i^K.rK (43) 

^This is only true within the one-flavor approximation considered here. As shown in ref. [50] things can 
change once flavor is considered. 
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Figure 5: The numerical solutions for a scenario with 10 decaying singlets that is in all 
other respects identical to the setup in in figure 4- The green/dashed line is the equilibrium 
abundance per comoving volume of one of the particle species, the red/solid line is its actual 
abundance as given by the Boltzmann equations, and the blue/dotted line is the absolute 
amount of the corresponding density of B — L. For the final value of the last quantity we 
find 1.2 ■ 10^"*^^ and therefore a slightly smaller value than in the corresponding one-particle 
scenario, which agrees with our considerations. 

The number density of each species of right handed neutrinos at this point is 



If all Ues singlets are to take part in the process, the final asymmetry will be the largest, 
when Ki = Mm{{Kj, I < j < Ues}) and Si = Max({e:j, 1 < j < ^eff}) for all i, which makes 
the different particle states identical within our treatment. However, we can also see that 
in this case the final baryon asymmetry does not change compared to the one particle case 
within our approximation, since we get a factor of Ues in the numerator and the denominator 
of eq.(45). A closer look at eq.(44) shows that the asymmetry will even be slightly reduced 
compared to the one particle case due to the factor z~^, which tends to smaller values. A 
numerical example that verifies our considerations is shown in figure 5, where we consider 
the same scenario as in section 2.2 with 10 identical neutrinos. 

3.3 The weak wash-out scenario in a radiation dominated universe 

As we have seen, in the strong wash-out regime the presence of many heavy neutrinos 
cannot lower the mass bound for heavy neutrinos compared to the standard case . This 
is due to the competing effects of the singlets on particle number and the wash-out rate. 
In the weak wash-out regime, this is different, as a result of the fact that an increase of 
the wash-out rate enlarges the final asymmetry. Therefore the increase of the number of 




(44) 



which leads to an approximate final B 



L asymmetry of 




(45) 
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decaying neutrinos and the increase of the wash-out rate, will not partially annihilate each 
other's effect. Instead both of them will amplify the final asymmetry. 

Let us make approximations similar to those in section 2.1. In concordance with equation 
(15) the abundance of each neutrino species with mass around the reheating temperature 
or lower is given by 

NnAz ^l)^ Hiz = l) ^ ■ ^^^^ 

The wash-out rate is the sum over the rates of all possible inverse decays and AL — 1- 
scatterings for which we can write 



2 

1=1 1=1 

in case of T > M,-. 

In the same manner in which we reached eq.(19), we now find 

"eff 



"eff / V \ 

Nb-l{z = 1) ^ ^£,iV^,(^^l)exp(^- ^^^'^^J 



(48) 



(J2^^K^)exp(^-J2K^) , (49) 



with the difference that now we are generally not allowed to make the linear approximation 
of the exponential function around zero. This is due to the fact the sum over the various 
Ki can be bigger than one, even though each Ki is much smaller. 

The late-time decay of the particles yields an asymmetry of YleiNfyf.{z ~ 1). Thus, com- 
bining inverse decay and decay periods, we are lead to the final asymmetry 

Nb-l(^ ^oo)^(J2 ^iKi) [exp (- J2 - l] ■ (50) 

In the case ^ i^Tj -C 1, this leads to 

Nb-l(z ^ oo) - ( J] SiK,) (X K,) , (51) 

while we can approximate 

Nb-l {z^(x)^-J2^iKi^ (52) 

in the case > 1^ 

Hence, it is evident that the presence of many right-handed neutrinos can easily increase 
the final B — L asymmetry in the weak wash-out regime. In the case of similar masses 
and couplings, the asymmetry increases proportional to for smaller and has a linear 
dependence on n^Q when J^Xj > 1. This is also illustrated by the numerical solutions 
shown in figure 6. 



^In the latter case, the phrase "weak wash-out" might be confusing, since the wash-out processes are in 
equilibrium around z « 1. However, around the decay time of the singlets, these processes will typically be 
frozen out. Therefore, we still refer to this scenario as a "weak wash-out scenario". 
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Figure 6: Here, we give two examples for leptogenesis with many singlets in the weak wash- 
out regime, where we used the same parameter values as in the setup considered in figure 
3. The red/solid line shows the abundance of particles of one singlet species in a comoving 
volume, whereas the green/dashed line illustrates the abundance of the corresponding equi- 
librium value. The blue/dotted line still represents the absolute amount of the corresponding 
B — L asymmetry. For the first picture we used n^jj = 10 and the absolute amount of the 
final asymmetry is found to be \NB~L{t — > oo)| = 1.0 ■ 10^"*^^, which is by a factor of 100 
larger than the corresponding one particle scenario, as suggested by our approximative ap- 
proach (cf. eq.(51)). In the second graph neff= 100 was used. Now, the wash-out processes 
are in equilibrium and the final asymmetry is found to be \NB-.L{t oo)| = 7.2 • 10~^°, 
which also roughly matches our estimation for this case (cf. eq.(52)). 

Of course these approximations are only valid, as long as neg does not get so big that the 
right-handed neutrinos get to dominate the energy density of the universe. In this case, the 
energy density of the singlets would become comparable to the background and they would 
produce a significant amount of entropy during their decay, hereby diluting the asymmetry. 
This possibility will be considered in the next section 

To get a feeling up to which point our estimations from this chapter are still valid, let 
us again assume that all of the decaying singlets have approximately identical masses and 
couplings {Ki = K If we postulate that the corresponding energy density never comes 

to dominate the universe, the strongest bound on the number of states will arise at the time 
of their decay. In this case, we can assume to be in a background dominated universe at 
early times, which implies that we can still use eq.(46) to determine the abundance of the 
various singlets. 

Once the states are populated, they will decay around the time 

T-^ = K-^H-\z = l), (53) 

where t=0 corresponds to z = 1. 
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Figure 7: Here, we show the energies stored in a comoving volume containing one photon, 
as calculated from the numerical solutions of eqs.(36) and (37) in a weak-wash out scenario 
(K = 0.01) with 900 different singlet states and parameter values as in the corresponding 
figures before. The energy due to radiation is illustrated by the green/dashed line, whereas 
the energy due to the singlets corresponds to the red/solid line. We can see that the energy 
densities become comparable around the time of decay, as suggested by eq.(59). 

By then, the radiation density will have decreased to 

TT^ - / aiz = 1) \ , , 

Pi = t;7T^?*^o TT^-FTTT (54) 



30^* \a{t = T-^), 
^ 35M^{2K~Y^ (55) 
^ 9M^K\ (56) 
where a is a comoving length scale. 

The energy stored in the right-handed neutrinos will also have decreased 

PN ^ neffK(^2^T3^M(2K-i)-^ (57) 

^ 0.1 ■ ries M^K^ , (58) 

with ^(3) ^ 1.2 [35]. 

The condition that Pn ^ Pi dX the time of decay, then translates into the relation 

90 

neff<^. (59) 

V K 



Also, our numerical analysis in figure 7 with 900 identical particles and K = 0.01 shows 
that the energy densities tend to become comparable around the time of decay if eq.(59) is 
fulfilled. 

3.4 Weak and strong wash-out: 

Before we consider the case of a singlet-dominated universe, let us stick to the background 
dominated universe a little longer and consider the case, where parts of our set of thermally 
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Figure 8: This is an example for a leptogenesis scenario with one strongly coupled right- 
handed neutrino (K = lO^j and ten weakly coupled ones (K = 10'^). All the other parame- 
ter values were chosen as in the examples before. The green/dashed line shows the thermal 
equilibrium abundance of each of the species in a comoving volume containing one photon. 
The corresponding actual abundances are illustrated by the red/solid lines, whereas the ab- 
solute amount of the B — L asymmetry in the same volume is shown by the blue/ dotted line 
and takes the final value 1.3 ■ 10^^'^. We see that the strongly coupled singlets closely follow 
the equilibrium value and wash out the B — L asymmetry produced during the population of 
the singlet states. At a later point the weakly coupled particles decay and the main part of 
the final asymmetry is produced. 

excited right-handed neutrinos have Ki <^ 1 and others Ki ^ 1. We will call the corre- 
sponding numbers of particle species ^ofr^w and riefr.s and further assume that all of these 
particles are within the same mass range. With the considerations from the previous pages, 
it is easy to guess, what we can expect from this scenario. Namely, while the presence 
of the weakly coupled particles will not have an effect on the final asymmetry produced 
by the strongly coupled ones, this cannot be said for the opposite case. The presence of 
one strongly coupled right-handed neutrino is enough to wash-out the asymmetry produced 
during the population of all the different states, including the weakly coupled ones. This 
implies that the asymmetry produced by the decays of the weakly coupled singlets is not 
canceled by an asymmetry of the same size but opposite sign, as it was the case in some 
of the previously considered scenarios. Hence, the final asymmetry in a scenario with one 
additional strongly coupled singlet can be substantially increased compared to a scenario 
with only weakly coupled particles of the same number. 

Combining our estimates from the previous sections, we can approximate the final asym- 
metry by 



where the sums in the two terms run over different particles, of course. In figure 8 we give a 
numerical example for one strongly coupled and ten weakly coupled particles, which nicely 
agrees with this estimate. 




(60) 
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3.5 The weak wash-out sceneirio in a matter dominated universe 



In equations (52) and (60), we have seen that we can increase the final baryon asymmetry of 
many leptogenesis scenarios by increasing the number of right-handed neutrinos. However, 
these estimations correspond to solutions of equations (36) and (37) and are therefore only 
valid in a universe which is dominated by background radiation and where the energy loss 
of this background due to the production of the singlet states can be neglected. 

However, if we have so many weakly interacting particles that the energy stored in this sector 
becomes comparable to the radiation energy of the background, their production rate starts 
to decrease even if they are far below their equilibrium value, since there is less and less 
energy left to populate these states. Additionally, these particles are massive and weakly 
coupled, which implies that the energy density of the radiative background will decrease 
faster than the energy density of the singlets. Thus, if the singlets live long enough they 
will come to dominate the universe prior to their decay. 

In these cases, the dynamics of the system are no longer governed by equations (36) and 
(37). Instead, we use the more complex system of equations 

dNj,^ 1 



dN, 



dt' Ho 



B-L 



Tri^D,n + rAL=l,n) (A^iV„ - N^) , (61) 
-no 

-JT E [^-{rD,n + rAL=l,n) (A^iV„ - KD + Tw^.n^^B-L] , (62) 



n= 



dt' Ho 

d{py) . 3.H{t') ^d{N, 



dt' '"^^ ' Ht, ^ dt' 

da H{t' 



-(P,«^)^-E^(^^»). (63) 

where t' = tHo is the time measured in units of Ho = LGGgl^"^ / Mp\. a is a comoving 
length scale with the volume containing one photon at early times (of course, the number 
of photons in this volume changes with time due to the production and decay of the singlets). 
Nn^ is the abundance of particles of the species Nn in a comoving volume of size a^, and 
Nb-l is the corresponding value for the B — L asymmetry. 

Further we have the relations 



<^") = ^^-^"HtM- '''' 

T = (^^^"Vt) ' , (67) 
where {E^) is the mean energy [32] of the Boltzmann distributed species Nn.^ 



^ Since the singlets are only weakly coupled to the SM lepton doublets and have no further interaction, 
the Boltzmann distribution can only be considered as an approximation. 
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However, it seems to make limited sense to consider arbitrarily large numbers of active 
singlet states ries, due to the fact that once the energy density stored in the right-handed 
neutrinos becomes comparable to the energy density of the radiative background, there will 
not be enough energy left to significantly increase the number density of the right-handed 
neutrinos. Therefore, an increase of ncfr beyond a certain limit does not imply a significant 
increase of populated singlet states. On the other hand, for very large values of rics the 
radiation energy might be transferred to the singlet sector so quickly that there will not 
be enough time for the wash-out processes to prevent the cancellation of the asymmetries 
during creation and decay of the singlets. 

Thus, we restrict our considerations to scenarios, where the number density of each singlet 
state per photon is still given by eq.(46), which can be quantified by the condition pNiT ~ 
Ml) < p-y{T f=:i Ml), which implies 

n,nKMt<^gM, (68) 

where we additionally assumed that all significantly abundant singlet states have similar 
masses around Mi. 

In this class of scenarios, the particle number density corresponding to each state at the 
respective time of decay will be 



[^-^—^^j ~ ^effi^ (^2 ^Mi^J [^K-^j ~ 0.1 • rieffX^Mf , (69) 

where we now need to determine the absolute number density, since the photon number in 
our original comoving volume will not be constant anymore. 

From this number density, we can get the energy density [p^ ~ MiUn), the approximate 
entropy density after the decay (s ~ Agl^'^p'^^'^/S [35]), and the number density of photons 



3 



n^f^2 ^^^J s ^OA-lQ-^ntKi . (70) 

' TT^ en 1 

Thus, our final value for the B — L asymmetry per photon (at early times) is 

f« -£i25n|ffi^3. (71) 

Th-y 

A numerical solution using eqs.(61)-(67) for this scenario is shown in fig. 9, which roughly 
matches our estimation for the final asymmetry. 



3.6 Lower bound for Mi 



As explained in the previous sections, the presence of many right-handed neutrinos will only 
help to increase the final baryon asymmetry (and hereby help to lower the lower bound on 
Ml) in the weak wash-out regime. Among the more closely considered scenarios, the one 
with the most decaying singlets will yield the largest final asymmetry. 
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Figure 9: An example for a leptogenesis scenario in a singlet dominated universe. For 
these two graphs, the parameters K = 10~^ and rieg = 3000 were used, with all other pa- 
rameters as in the previously considered cases. The left-hand side shows the equilibrium 
abundance of the singlets (green/ dashed) per photon, their actual abundance (red/solid) and 
the absolute amount of the B — L asymmetry (blue/dotted) per photon which is found to be 
8.4- 10~^ at the end of the calculation, which roughly agrees with our estimation in eq.(71). 
The time t was chosen such that t = corresponds to the singular point in the previously 
radiation dominated universe. The graph on the right-hand side, shows the radiation energy 
in a comoving volume that contains one photon at early times (green/ dotted) and the corre- 
sponding energy stored in the singlets (red/solid) in GeV. We can see that the singlets get 
to dominate the universe before their decay. 

Therefore combining eq.(68) with eqs.(5),(8), and (71), we find 




With K 10^^, we find an approximate lower bound around of ~ Mi > lO'^GeV. An 
increase of K by an order of magnitude or so, would slightly increase the result of this 
formula. However, the numerics show that we would not be allowed to neglect the wash- 
out processes during the decay anymore and that the final asymmetry already tends to be 
smaller than in the K = 10~^ case. 

Numerically we were able to push the lower bound a little further as we were able reproduce 
the observed BAU with K = 10~^, Ues = and 

T^ = Mi = 6.5 ■ lO^GeV. (73) 

We see that this relaxes the lower bound in eq.(32) for standard thermal leptogenesis by 
approximately one and a half orders of magnitude. Moreover, the relaxation of the Davidson- 
Ibarra bound (as discussed in sec. 3.1) as well as a complete parameter scan might motivate 
a further relaxation of this bound. 

3.7 Upper mass bound for 

The analysis in section 2.3 that lead to the upper bound for in the standard case is 
not valid in the presence of many singlets. The reason for this is that we cannot draw the 
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conclusion that we are in the strong wash-out regime from the fact that we have degenerate 
hght neutrinos as we did in the standard case. All we need now is that three singlets are 
in the strong wash-out regime. However, if we have many more singlets that are weakly 
coupled, they can determine the size of the final baryon asymmetry and not the strongly 
coupled ones that might determine the size of the light neutrino masses. 



Thus, even if we imply the Davidson-Ibarra bound, we find that within our rough approxi- 
mations the bound on the absolute neutrino mass scale becomes less stringent than bounds 
from other areas (e.g. [51-53]). In scenarios that evade the Davidson- Ibarra bound, the 
upper mass bound on the light neutrino masses can be relaxed even further. Additionally, 
the bound in eq.(74) seems rather conservative from the point that the AL = 2 scatter- 
ings do not need to be frozen out around z — 1 in the weak wash-out regime but at lower 
temperatures. 

3.8 Further Issues 

AL = 2 Scatterings: 

So far we have neglected the possibility of AL = 2 scatterings. Our estimate for the 
corresponding bound, is the same as the one wc used in the standard thermal leptogenesis 
case. As the Davidson-Ibarra bound, it does not get modified by the presence of many 
singlets, since the increase of the number of relevant diagrams is compensated by the increase 
of diagrams used for the see-saw mechanism. Again, we estimate the same reaction rate for 
the corresponding processes as in sec.2.5 and therefore also use eq.(33) as a criterion for the 
validity of our considerations. 

Quantum Effects: 

As already mentioned in our review of standard thermal leptogenesis in the previous section, 
there arc also questions about the possible contributions of quantum effects. However, a 
detailed analysis of these effects in the presence of many singlets goes beyond the scope of 
this paper. 

Initial Conditions: 

While the initial conditions for the singlets abundances are not typically important in the 
strong wash-out regime, they are a crucial point in the weak wash-out regime, in which 
we are more interested. So far, we have always made the conservative assumption of zero 
initial abundance of the singlets. Without being quantitative let us make some comments 
on further possibilities. 

If we assume to have a thermal initial abundance of the singlets (democratic reheating), we 
should also be able to drop the corresponding lower bound for the reheating temperature 
from refs. [7,8] significantly, as a large number of singlet states would automatically imply a 

universe dominated by singlets whose decay produces a lepton asymmetry, cf. [8]. However, 
compared to a universe dominated to by one singlet species we would probably not be 
able to lower the bound for the reheating temperature, since the only difference (in case 



Combining eqs. (29), (30), (8), and (71) we find 




(74) 
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of identical states) would be an increase of the wash-out rate, which might even slightly 
decrease the final asymmetry. Thus, while the setup also seems interesting for democratic 
reheating scenarios it seems less interesting for non-thermal baryogenesis scenarios, e.g. [54]. 



4 An Extra-Dimensional Example 

In this section we present a possible realization of the previously discussed leptogenesis 
scenario within a concrete model. Therefore, we adopt the frequently used idea that gauge 
singlets can propagate in an additional dimension, while all the standard model particles are 
confined to a 3-brane [55]. If there are no further interactions, one can formally integrate 
over this extra dimension and is left with a regular four dimensional theory with a Kaluza- 
Klein tower of right-handed neutrinos. For a broad parameter space, this can lead to a 
quasi degeneracy of many of the Kaluza-Klein modes, which is a key ingredient needed for 
leptogenesis with many singlets, as we illustrated in the previous section. 

In this context, let us mention ref. [56], which also introduced an extra- dimensional model 
for a leptogenesis scenario that made use of the Kaluza-Klein tower of extra-dimensional 
neutrinos. However, the model is quite different from our proposal in many aspects. This 
is partially due to the different nature of the 5-d Majorana mass term, which is scalar- 
hke in our model and vector-like in ref. [56]. Another difference is that ref. [56] needs 
to violate the orbifold symmetry to get a non-zero CP violation (since it only uses one 
generation of singlets), which is not necessary in our case. Also, we do not make use of 
any resonance effects that increase CP violation. Finally, since ref. [56] assumes much 
larger extra-dimensions, its leptogenesis scenario takes place after the electro-weak phase 
transition, which makes the quantitative treatment more difficult. The model from ref. 
[56] has been considered further in ref. [57], which also suggests the consideration of three 
generations of singlets in the bulk to increase the final asymmetry. 

Refs. [58,59], e.g., also consider extra-dimensional leptogenesis scenarios. While ref. [59] 
state that many decaying Kaluza-Klein states can increase the final asymmetry of their 
scenario, ref. [58] finds that a possibly decaying Kaluza-Klein tower of particles would not 
increase the final asymmetry in the respective model. 

4.1 The model 

Let us be more quantitative and use the treatment from refs. [47,55,60,61], where the 
groundwork of this subsection has been developed. 

As already mentioned, we use a five dimensional model, where the standard model particles 
are confined to the regular 3-brane. Using the notation of ref. [62] for this section, the 
corresponding part of the action is given by 



Here, y is the additional spatial degree of freedom, which we assume to be compactified on 
an orbifold Sijlji of radius R. 




(75) 
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Additionally, wc introduce three singlet fields that can propagate in the extra-dimension 
and which transform as P^'^i = 75 under the parity transformation : y ^ —y. This 
allows us to write down the gauge and orbifold-symmetry invariant action 



'bulk 



/ 



d!^x dy 



*,i7°9«*,--(M,,*f*,+h.c.) 



where is the corresponding charge conjugate field of defined by 



e 
e 



(76) 



(77) 



with e = ia2. We can always choose the such that — diag(Mi, M2, Ms)ij with Mj e R 
and < Ml < M2 < M3, which we will therefore assume in the following. 

Further, we can introduce brane-bulk couphngs that also respect the SM gauge invariance 
and the orbifold symmetry 



^1 



brane— bulk 



d'^x dy 



9ij 



■.■^iPLH {{0,0) Jjf + h.c. 



(78) 



where the Qij are dimensionless Yukawa couplings and M5 is a mass scale introduced for 
dimensional reasons. 

The transformation properties of the allow us to decompose them into modes 

[ ^yZi^) + £ ^os{ny/R) *g(x) \ 

n=l 

ZMny/R) 41 (a:) 



'^i{x,y) = 



\ 



(79) 



n=l 



with the respective right- and left-handed Weyl spinors ^^^^ .(x). 

We can now perform the integration over y in the action and are left with a four dimensional 
theory with infinitely many particles due to the Kaluza-Klein tower of the ^j. After the 
Higgs field develops a VEV {(f)) we find an infinite Majorana mass matrix for all the left- 
handed neutrino fields {u, ^f^, "^ff, . . . ) 



/ -ml/V2 
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-M 
-n/R 



-m 
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-n/R 
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\ 



(80) 



with = -e[¥^l]* and mD,ij = gij{(P) /^/Mf,. 
If we further diagonalize the sub-matrices 

-M -n/R 
-n/R M 



(81) 
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by a change of basis {u, ^f^, ^'^ , . . . ) - ^J?, ■ ■ ■ , ^f^J, ^J?, . . . ), the Ma- 
jorana mass matrix becomes 



/ 



-mn,+ 



M 














^jM^ + v?/R^ 



\ 



(82) 



where 



and 
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(83) 
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We see that the Kaluza-Klein tower consists of pairs of fields with degenerate masses. Ad- 
ditionally, in case of M;^R ^ 1 a large number of pairs pairs will be quasi-degenerate with 
respect to their masses. Also, we note that the relative strengths and phases with which a 
mode couples to the SM neutrinos is the same for all the modes. Still, the couplings of the 
neutrinos to the different Weyl fields can differ by a constant factor and a possible overall 
phase. 

Let us take a closer look at the Yukawa couplings of the various modes ^j, namely 



(^+,n)jj — 



(0) 



ITTlr 



(0) 
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9ij 
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^Mf + nVi?2 ' 



y/2nRM5 
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(86) 



(87) 



^The author would like to thank Naoyuki Haba for related, private notes on the five dimensional see-saw 
(with vector- like masses) [63] that inspired this treatment. 
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for n < 1. We see that for n ^ Mji? the couphng strengths of the ^j^l and modes 
become similar, whereas for n -C MiR we find 




n 



(Y_\.. ^ ^ (89) 

from which we see that the couphngs of the lower modes of ^^^^^^ can be strongly suppressed 

with respect to the lower modes '^^^\ if the corresponding Majorana mass scale Mj is much 
larger than the compactification scale R~^. 

4.2 CP asymmetry 

In the context of leptogenesis, it is of course important to determine the CP asymmetry 
produced by decaying modes of ^j. With our considerations from previous chapters, this is 
now straightforward/ 

Applying eq.(40) we find 

with Mj_„ = Mf + rfi jW-. After further transformations we get 

Thus, each mode features two particle species, which yield exactly opposite CP asymmetries. 
Further the CP asymmetries only depend on the five dimensional wave number through the 
proportionality factor Mj_„. 

Let us also stress that in the case M3 > M2 » Mi » in which we have many 

quasi-degenerate neutrinos, we do not get any resonance or higher-order effects as discussed 
earlier, since all the quasi-degenerate particles with masses around M\ have the same Yukawa 
couplings (except for a purely real or imaginary overall factor). 

Finally, we note that all singlets of one generation have the same Yukawa couplings (up to 
an overall factor) . This leads to the fact that the effective light neutrino mass matrix can be 
written as the product of 3 x 3 matrices instead of n^r x 3. E.g. for MiR ^ 1, the neutrino 
mass matrix in our model is given by [47,55,60,61] (see also [62,64]) 



''Since there are cancellations between the various contributions to the see-saw in this model, the per- 

turbative approach using Feynman diagrams might lose its validity. However, if the cut-off for our extra- 
dimensional theory is around the scale M5, this problems seems to be absent. 
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Yet, this yields the relation 



In this case the Casas-Ibarra parametrization [49] of the Ynkawa couplings can be applied 
(with a possible additional overall factor and the heavy masses being replaced by 2/{t{R)). 
After some easy algebra with eq.(40) one then finds that the more stringent Davidson-Ibarra 
bound given in eq.(5) is still valid in this scenario. 

4.3 Leptogenesis 

In addition to the masses of the various modes and the CP asymmetries they yield, the 
decay rates of the modes of the first generation are also needed to determine the final baryon 
asymmetry. Similarly to the previous chapters it is useful to characterize this quantity in 
terms of the parameter 

_ r^,.,^(T ^ 0) [(y^,n)(y^,n)l,M,,„/(87r) / M,\ 

H{T^M,) - IMgMl/M^, " I M^, J ' ^^^^ 

for n > 1 with a corresponding definition for Kq 

Considering eq.(93) and the four-dimensional effective neutrino mass (cf. eq(ll)), one can 
define the five-dimensional effective neutrino mass as 

m^M = (m£,m d) ii ^ ■ (95) 
m* ttMiR 

which implies that the parameter Kq is pushed to smaller values with respect to the corre- 
sponding four dimensional scenario. Hence, while the size of the measured neutrino mass 
differences might make K ^ 1 more plausible in four dimensions, the weak wash-out regime 
can be more natural in this scenario. 

Once we have fixed the values Mi, R, Tr, and Kq we can already use eq.(52) and (71) 
to approximate the final B — L asymmetry of this leptogenesis setup. The key feature of 
this model is the fact that modes that yield opposite CP asymmetries couple with different 
strengths to the lepton doublets. Being in the weak wash-out regime, their contributions to 
the final B — L asymmetry differ in size, which prevents them from canceling each other. If 
the reheating temperature equals the mass of the lightest of the heavy neutrinos (Tr = Mi), 
there are of the order of MiR quasi-degenerate neutrinos which will reach significant particle 
abundances. Thus, if we assume maximal CP violation, a rough estimate for the final B — L 
asymmetry is 

Nb-l ^ 10^^ i ( -4^) {M,R) K, , (97) 
Vo.OSeV/ VlO^GeVy ^ i ; O' v ; 

or 

^-'-^•^■^°"(ai5v)(To^)(^'«)'*^ (^«) 

depending on the fact if the singlets come to dominate the universe or not, and therefore 
on the size of M^R (cf. eq.(59)). 
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Figure 10: The two graphs show the numerical solution for MiR pairs of particles with 
parameter values as described in the text. In the graph on the left, the green/dashed line 
represents the equilibrium abundance of each of the particle species per photon. The red/solid 
line shows the corresponding value for a "+"-mode, while the behavior of the corresponding 
"-"-mode is illustrated by the yellow /dashed- dotted line. The absolute amount of the B — L 
asymmetry is again represented by the blue/ dotted line and reaches a final value 0/ 7.8 -10^^. 
In the second graph the behavior of the respective energies in a comoving volume is shown, 
where the red/solid line corresponds to the singlets and the green/dotted line to radiation. 
The comoving volume contains one photon at early times. 

Using eq.(8) we can see tliat values sucli as Mi = 2-10® GeV, = 10^ GeV, and Kq = 10"^ 
already seem to yield a baryon asymmetry of the right order of magnitude. 

Let us be more precise. If equation (52) and therefore also eq.(97) is vahd, we see that the 
contribution of a pair of particles \ to the final baryon asymmetry is independent of their 

actual masses, since K^+l - i^i"! = 2Mi/Mi,„, which means the Mi „ in the denominator 
will cancel with the Mi„ in si. Thus, all pairs with masses low enough that there will 
be a significant abundance of them approximately contribute the same amount to the final 
asymmetry. 

If equation (71) and therefore also eq.(98) is valid, we need to consider {K^\)^/^ — {K^\)^/'^, 
which means that pairs with larger masses will contribute slightly less, even if they get 
produced in a significant abundance. This means, instead of treating all MiR particle pairs 
as different, we can work with a system of completely similar particle pairs of mass Mi^„ 
to get an approximate lower bound for Mi and the reheating temperature. This has been 
done in figure 10 numerically, where we were able to produce a. B — L asymmetry that was 
even slightly larger than the observed value with = 10^ GeV, Kq = 10~^, and 

T, = Ml = 2 • lO^GeV, (99) 

while we also assumed maximal CP asymmetry. 

Of course, lower values would be possible if we would consider also < lO^GeV. However, 
due to our arguments in the next subsection this does not seem to be motivated. 
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4.4 Constraints From Extra-Dimensional Gravitons 



In addition to the open points listed in section 3.8, this extra-dimensional scenario has 
to consider the effects of gravitons that are likely to propagate in the extra dimension. 
In particular, we need to ensure that excitations of gravitational modes do not spoil the 
leptogenesis scenario by producing too much entropy, especially during BBN. Therefore, 
we consider the possible production and decay of gravitational modes, mainly using the 
analysis from ref . [65] . 

Since the coupling constant of a 5-d graviton to the standard model particles is proportional 
to l/Mg^^, a single mode coupling is proportional to l/(M|i?)^/^ (cf. eq.(79)). By dimen- 
sional analysis the decay rate of a single Fourier mode of mass Mq,. in its rest frame (from 
the four dimensional point of view) is therefore 

r„,c,-.M„^ = ».o^, (100) 

where ride is the number of decay channels. 

If we define the parameter Kq,. for a graviton mode corresponding to eq.(lO), we get 

The typical decay time of these modes is Tq] = Kq]H{T — Mq^)~^ (cf. eq.(53)), which 
corresponds to the temperature 



■Gr 



^Gr ( f ) = Mo.{2K-^])-i (102) 




> 



10 TeV 



3 



(103) 

(104) 

• 0.3 MeV. (105) 



where we estimated Udc ^ g* in eq.(104). 

Since the lowest mass of a gravitational mode (other than the zero mode) approximately 
equals the compactification scale, we find the requirement that the compactification scale 
approximately needs to be of the order 50 TeV or larger. Otherwise, the gravitational modes 
might decay late and change the predictions of big bang nucleosynthesis. 

Even if the gravitons decay early enough to avoid bounds from BBN, it is also important that 

they do not get to dominate the energy density of the universe before that time. Otherwise 
the entropy production during their decay could invalidate our estimations for the size of 
the produced baryon asymmetry. 
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To consider this possible issue, we can use eq.(46) to roughly estimate the abundance of a 
gravitational mode in a comoving volume, which contains one photon. Since the number 
of modes with mass Mqj. is approximately Mqj.R we can estimate the energy per comoving 
volume stored in gravitons of mass Mqj. at the time of their decay by 

EGr ^ (MarR) Kor Mq. ^ , (106) 

whereas the radiation energy stored in the same volume at the same time is approximately 
g,TGj2 (cf. eq.(104)). 

The ratio of the two energies is therefore approximately 



which indicates that even for a compactification scale as low as 10 TeV, the reheating 
temperature can be as high as 10^° GeV without a significant entropy production due to 
gravitons. 

Thus, within the limits of our estimations the lower limit on the reheating temperature 
for the leptogenesis scenario considered in this section does not get additional bounds from 
the consideration of gravitational modes other than a compactification scale in the 50 TeV 
range or larger. 



5 Summary and Outlook 



Considering leptogenesis scenarios with many singlet neutrinos, we have found that it is 
possible to lower the reheating temperature of standard thermal leptogenesis by approxi- 
mately one and a half orders of magnitude. This effect can be achieved in the weak wash-out 
regime, while the strong wash-out scenarios become slightly more restrictive compared to 
the standard case. We also argued that the weak wash-out scenario can become more nat- 
ural in theories with many singlets. Approximative formulae have been given that help 
to illustrate the effect of the additional particles. Additionally, we have checked our ap- 
proximations by comparing them to results from numerical solutions of the corresponding 
Boltzmann equations. 

We further argued that the bound on the absolute neutrino mass scale from standard thermal 
leptogenesis can get weakened in the considered scenarios and might fall behind bounds from 
other areas of particle and astro-particle physics. 

As a specific realization of a leptogenesis with many singlets, we considered an extra- 
dimensional scenario, where the Kaluza-Klein excitations of right-handed singlets serve as 
right-handed neutrinos in the weak wash-out regime and provide all necessary ingredients 
for leptogenesis below the usual bounds. We also considered the excitations of gravitational 
bounds for this model. Within our approximations we find that these bounds still allow for 
a significant reduction of the necessary reheating temperature. 
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Neither for the general case of many neutrinos, nor for the exphcit extra-dimensional model 
we have explicitly performed parameter scans. This means, there might be a parameter 
range for which the bound for the reheating temperature can even be lowered further. 

The fact that the Davidson-Ibarra bound docs not generally hold in the case of many 
singlets might also leave room for a further relaxation of the lower bound for the reheating 
temperature. This is especially interesting in the context of our arguments for the lift of 
the upper mass bound for the light neutrinos masses. 

An interesting topic that goes beyond the scope of this paper is the inclusion of flavor 
dependence in this setup, especially, since it may be that the interference of tree-level and 
the self-energy correction diagrams for the CP asymmetry gives important additional effects, 
if the particles in the loop propagate in the opposite direction compared to fig. 2. 

6 Acknowledgments 

I would especially hke to thank Yosef Nir for valuable comments on the first version of this 
paper, in particular in the context of the Davidson-Ibarra bound with many neutrinos. 

Further, I would like to thank Michael Ratz for numerous fruitful discussions that had a 
large impact on the quality of the present work, as well as for valuable comments on the 
manuscript. I would also like to thank Naoyuki Haba for introducing me to the idea of extra- 
dimensional neutrinos, discussions on the 5d see-saw and leptogenesis at early stages, as well 
as for private notes on the 5d see-saw. Moreover, my thanks go to Michael Pliimacher and 
Oleg Lebedev for useful discussions. Thanks for financial support go to the "Graduiertenkol- 
leg 1054" and the "SFB-transregio 27" of the "Deutsche Forschungsgemeinschaft" as well 
as the cluster of excellence "Origin and Structure of the Universe" . 



30 



A Reaction Rates 



For our numerical solutions we used the reaction rates as mainly listed in ref. [9]. 
The decay rate is given by [32] 

where the Ki{z) are modified Bessel functions of the second kind. 

The AL = 1 scattering rates for the s- and t-channel processes are given by [5] 

Ss ^ ^-^^^2.^Uz) 
Hz 12-"^ ' 

St ^ ^^ = 2-2-%,(z) 
Hz 12-' ^ ' 



(108) 
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where we had to introduce a Higgs mass (o/j = (m/j/Mi)^) to regularize the infrared diver- 
gence of the t-channel diagrams but treated the Higgs boson and all other particles (except 
for the right-handed neutrinos) as massless everywhere else . 

To cut down the time needed for the numerical calculations we only calculated fs/t{z) at a 
fixed number of points (20-30) and then linearly extrapolated the functions. 

Further, the wash-out rate is given by 
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where A^tvi is the actual abundance of the right-handed neutrinos in a comoving volume 
containing one photon, iV^ = 82:^/^2 (2^) /8 is the corresponding value for an equilibrium 
distribution, and A^^^'* = 3/4 is the corresponding value for massless leptons. 

The generalization to the case of many decaying the neutrinos is simple (cf. eqs. (36) and 
(37)), since all one has to do is to do the interchange Mi M„ and K — > for each singlet 
state. ^ However, one has to remember that the factor in the denominator Hz does obviously 
not change for each state. To find the correct rates one has to make the transformations 

and in the same way for the S and W terms (also ah — > a/j(Mi/M„)^). 

With these considerations it is also simple to get the rates for the system described by 
eqs. (61) to (67). Namely, 

^^^dMlllIl ^ D^{M^,T) ^ (119) 
Hq Ml 

and accordingly for S and W. 

B The Davidson-Ibarra bound revisited 

In this part of the Appendix we generalize the treatment from ref. [16] to find an upper 
bound for the CP-asymmetry for a setup with 11^ heavy singlets and three light neutrinos. In 
this case, one can generalize the Casas-Ibarra parametrization [49] of the Yukawa couplings 
to 

9 = ^D^RD^U^ , (120) 

where -Dy^g and are diagonal matrices that contain the masses of the heavy and light 
neutrinos, respectively. Here, we work in the basis where the heavy mass matrix is diagonal 
and U is the mixing matrix of the light neutrinos, i? is a complex x3 matrix that fulfills 
the condition 

R^R = 13X3- (121) 

Let us use a mathematical trick and introduce n^v — 3 additional light neutrinos with zero 
couplings. We can easily do this, since these particles do not show up in any particle physics 
experiment. Now, the Yukawa couplings can be parametrized by 

9 - ^.D^R'D'^U'^ , (122) 

with = diag(mi, 1712, rris, 0, . . . , 0) and 

\^ U l(n^-3)x(nAr-3) / 



^It is important to note that, aside from the change due to the Yukawa couphngs, Kn acquires a factor 
of Mn/Mi and not of the inverse, since only the Mi factor in the decay rate gets modified while Hq does, 
of course, not change. 
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R' is now an arbitrary complex orthogonal un x nN matrix with 

r>/T T>l -n 

As in the 3x3 case [16] eq.(40) can now be transformed to 

3 M, T.jm]l^{R%) 

with uik = for k > 3. 
Using Rfj = 1 we then get 

3 Mi Max({m2}) - Min({m2}) _ 3 
~ 167r (0)2 Max({mj}) + Min({mj}) ~ IGtt (0)2 

where rris is essentially a free parameter. 
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